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POTENTIAL-PLANE ANALOG SIMPLIFIES 
DESIGN OF NUCLEAR REACTOR CONTROL 


The potential-plane analog is based on the principle of using a voltage distribu- 
tion on a uniform resistive sheet to represent a function of a complex variable. The 
equipotential lines may represent the log of the magnitude or the angle of a function 
of a complex variable. If two resistive sheets are used and the electrodes on one 
are orthogonal to the electrodes on the other, both the log magnitude and angle func- 
tions can be obtained simultaneously. 


Although the principle of the potential-plane analog has been known for some 
time, it has seen little use in the past due to some severe limitations that made it 
difficult to construct a practical machine. An extension of potential-plane analog 
- theory by M. L. Morgan has eliminated the previous restrictions with the result 
that a practical potential-plane analog is now available. 


The theory and design of a practical potential-plane analog has been described 
previously (1,2,3) and the purpose of this paper is to show how a potential-plane 
analog can be used to simplify the analysis and synthesis of control systems, in 
particular, a nuclear reactor control system. In the following discussion, a nuclear 
reactor is described by a transfer function. This transfer function is then included 
as a component of the control system to be analysed with an ESIAC* potential-plane 
analog. Although simple transfer functions are used in this example to keep from 
drowning the analysis in a sea of mathematics, this should not be taken as a limita- 
tion on the ESIAC potential-plane analog, since it is capable of handling approxi- 
mately twenty poles and zeros. 


REACTOR PHYSICS 


A nuclear reactor is a system usually consisting of a moderator, fissionable 
fuel material; a means of removing heat, and a geometric structure suitable for 
maintaining a chain reaction. Figure 1 shows the basic components of one type of 
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power reactor. The fuel is located in the center section of the reactor and consists 
of uranium 235, in solid form. The moderator is intermingled with the fuel and is 
used to slow down the neutrons resulting from the fission process to where they 
possess the desired energy range. For this example, the moderator material is 
graphite and the neutrons are slowed down to the thermal energy range (0.026ev). 


The control rods are made of a neutron absorbing material such as cadmium, 
boron, or hafnium and they control the chain reaction by absorbing a percentage of 
neutrons produced in the fission reaction. In each fission reaction, the collision of 
one neutron with a U“”’” nucleus will produce approximately Z.5 neutrons that are 
each capable of causing other fission reactions. However, a percentage of these 
neutrons are lost for various reasons and the control rods are used to absorb the 
excess so that the average number of neutrons available to produce another fission 
is approximately one. 


The energy produced in the reactor core in the form of heat is transferred by 
the cooling system to an external plant which usually converts it to electrical energy. 
Outside of the reactor core is a reflector that helps to maintain the chain reaction 
by reflecting the neutrons back into the core that otherwise might be radiated and 
lost. A biographical shield surrounds the reactor and attenuates harmful radiation 
such as gamma rays and neutrons. The last basic component is a neutron detector 
which is used to monitor the level of neutron concentration in the core. A common 
type of neutron detector is an ionization chamber which gives an output current pro- 
portional to the number of fissions occuring per second. The power generated by 
the reactor is also proportional to the number of fissions per second. The relation- 
ship is: 


3.10'"° fissions/sec = | watt (1) 


Therefore, the neutron detector indicates the power level at which the reactor is 
operating. 


For any specific reactor, an effective multiplication factor k can be defined as 
the ratio of the number of neutrons in any one generation to the number of corre- 
sponding neutrons of the immediately receeding generation. Then, if kis equal to, 
or slightly greater than unity, a chain reaction can take place. Ifkis less than 
unity, a chain reaction cannot be self-regenerative, and will ultimately die out. In 
most reactor operations, kwill be very close to unity so another quantity, reactivity, 
is usually used. Reactivity Ak is defined as: 


Ak = a ~k—!l (2) 


Thus, reactivity is the amount by which the multiplication factor of a specific reac- 
tor differs from unity. 


Another important reactor parameter is the neutron lifetime T which is defined 
as the average effective time between successive neutron generations in a specific 
reactor. Since the number of neutrons changes from one generation to the next by 
Ak, and assuming that there are initially n neutrons per cubic centimeter, the rate 
of increase in each generation isnAk. The rate of increase with respect to time is 
then: 
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dn. Ak 3 
at fT n neutrons /cm”/sec (3) 


Integrating this equation gives: 


Nn =N€ Satis neutrons /cm° (4) 


where n, is the initial concentration of neutrons, and. is the concentration after 
time t. Equation (4) shows that the neutron concentration rises exponentially with 
time if Ak is positive which is the case for k greater than unity. Notice that the 
reactor period or time constant {_L\ depends upon the effective lifetime T and the 
reciprocal of the reactivity Ak. 


Fortunately for the control system designer, the exponential response given in 
Equation (4) is modified by the fact that not all the neutrons created in the fission 
process are given off instantaneously. Approximately 0.75 percent of the neutrons 
are delayed, and their average lifetime is 10 seconds. Actually, these delayed 
neutrons are given off in six distinct groups at different times and in different quan- 
tities. The symbol 8 is used to denote the total fraction of delayed neutrons (0.0075) 
and #j is the fraction of delayed neutrons in the i'th group. Similarly, the decay 
constant, which is the inverse of the avereecc Le eiwer is denoted by Aj with the 
average decay constant \ equal to 0.1 sec-1+. The presence of the delayed neutrons 
imparts a desirable "inertia" effect to the exponential characteristic of the reactor. 
Consequently, it is practical to control a reactor automatically compared to the 
very difficult problem it would be without the delayed neutrons. 


The instantaneous state of the reactor is described in terms of the value of the 
effective multiplication factor k. If k= 1, the reactor is said to be critical k <| 
subcritical, and k >!supercritical. If k= 1.0075, the reactor is said to be prompt 
critical. This means that the reactor can sustain a chain reaction without the use 
of the delayed neutrons. Generally, most control systems prevent k from becom- 
ing greater than 1.0075 because the reactor characteristic becomes a rapid expon- 
ential when k is greater than 1.0075. 


REACTOR TRANSFER FUNCTION 


In the eyes of the control engineer, a reactor is a black box whose character- 
istics are specified by a transfer function. The reactor transfer function can then 
be incorporated with the transfer functions of the other system components to get 
the over-all system transfer function. 


Since the behavior of the over-all system is specified by the system transfer 
function, the control engineer can then determine whether or not the system works 
properly, and if not, he can modify it so that it does. 


It has been shown that a reactor can be described in the time domain by the set 
of differential equations shown below: 


dn _ Ak—8B SA 
“dt. - T n+ 2 ri (5) 
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where fj; is the number of delayed neutrons per cubic centimeter emitted in group i 
and the other symbols are the same as defined previously. For control purposes, 
the reactor is best described by a function in the frequency domain and this function 
Ak may be found from the above equations with the aid of the Laplace Transform 
a shown in the appendix. The output quantity of the reactor transfer function is a 
change in neutron concentration An and the input quantity is a change in reactivity 
Ak , which is due to a change in control rod position. The reactor transfer function 


is derived in the appendix and is: 


An(s) _ l 


Ak(s) = T - B; 
sf +) T(s +i) 


(7) 


and by expanding the denominator, 


ante) ny UL +N) 


ee, ee Oe (8) 
Ak(s) U s TI (s +r;) 


where s=o0+jw is the Laplace Transform variable andw=27f. Using the values of )j 
and B;, given in reference (5), and a lifetime T of (T=lO0;gconas)» the poles r; were found 
by factoring the sixth degree polynomial that results from expanding the denomina~ 
tor of Equation (7). The result is a reactor transfer function of: 


An(s) _ no (8 +14)(s+1.61)(s + 0.456) (s+ 0.15!) (s + 0.0315) (s + 0.0124) 


Ak(s) 10% (77)(s+13.4)(6+1.43)(s+0.34)(s+0.08)(s+0.015) 


where n, is the steady-state neutron concentration and An is a change from the 
steady-state value. 


Figures 2 and 3 show the Bode plots (magnitude and Angle versus log f) of the 
above transfer function as plotted with an ESIAC potential-plane analog. They are 
normalized with respect to A with (T =lO;éconas)* 


Although the reactor transfer function shown in Equation (9) is not difficult to 
handle with an ESIAC potential plane analog, a simplified approximation suitable for 
most oontrol engineering is: 


An _ No (s +) 


e (10) 
Ak \o~* s(s+r) 


The values of A, \ and r can be chosen to give Bode plots that resemble those of 
Equation (9) as accurately as possible. 
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The transfer function: 


Ss 
An _ . 7l03)(s+0.07) _—‘|4 (1+ 533) 
Bk * " sle435)  g/ya 5) 7 
I ae 


was chosen for this paper and the Bode plots of Equations (9) and (1!) are compared in 
Figures 2 and 3. 


The elementary reactor transfer function discussed above ignores two funda- 
mental processes that are also needed to describe reactor operation. These proc~ 
esses produce a poisoning effect and a temperature coefficient effect. The poison- 
ing effect will not be considered in this paper because it is a very low frequency 
phenomena and is usually most important only during start-up of the reactor after 
a period of low power operation. 


The temperature coefficient effect is primarily due to a decrease in density of 
the moderator and reflector material as the temperature within the reactor is 
raised. The change in density of the moderator and reflector usually causes a 
decrease in reactivity as temperature increases and thus the temperature coefficient 
is said to be negative. The temperature coefficient is one of the most important 
reactor parameters and it is usually negative for a graphite-moderated, Us35 
thermal reactor. 


The temperature coefficient can be included in the over-all reactor transfer 
function by adding a feedback loop to the elementary reactor transfer function as 
shown in Figure 4. 


ELEMENTARY 
REACTOR 


TEMPERATURE 
COEFFICIENT 


FIGURE 4. BLOCK DIAGRAM OF REACTOR WITH 
TEMPERATURE COEFFICIENT 


The reactor transfer function including temperature coefficient can be found by 
using the feedback equation: 


OUTPUT _ G(s) 


SOUTEUT, 6B (12) 
INPUT 1+ G(s)H(s) 


and is: 
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aes 

Oe I4 no (1+ 533) 

Ak — s s 
s(\+ 35) + Kre 14n, ('+ 57) 


Intuitively, it would seem that the negative feedback due to the temperature 
coefficient would tend to stabilize the reactor. The effect of the reactor negative 
temperature coefficient on an automatic control system including a reactor is 
shown in the following example. 


(13) 


ILLUSTRATIVE EXAMPLE 


A block diagram of a nuclear reactor with its associated system components is 
shown in Figure 5. 
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FIGURE 5. NUCLEAR REACTOR AND CONTROL SYSTEM 


The output of the system is power which is proportional to the neutron concentra- 
tionn. The neutron concentration consists of a steady-state level n, and a variation 
An. The input to the reactor is reactivity Ak which is determined by the position of 
the control rods D. The position of the control rods is determined by the control 
rod servo which is driven by an error amplifier. The error voltage V, is the output 
of the comparator circuit and is described by: 


ae (14) 
e ~ V; 

This function gives an error voltage that is inversely proportional to steady-state 

power level n,, since Vj =kgn;, and nj is approximately equal ton,.. This inverse 

relationship is deliberately used to cancel out the effect of the n,in the reactor 

transfer function. Unfortunately, the amount of cancellation depends upon the 

temperature coefficient K;,as shown later, 


We have said previously that the reactor temperature coefficient is likely to be 
negative for a graphite-moderated thermal reactor. To begin the analysis of the 
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control system, let us assume that the temperature coefficient may vary from zero 
to some negative value and that the system will be least stable for the zero value. 
We can find the stability characteristics of the system for this worst condition and 
then find the change in the system due to the temperature coefficient becoming neg~ 
ative. 


The complete system transfer function with K;, =O is derived in the appendix 
and is: 


n _ 1+ G(s) 


Ti st-+G(s) H(s) = 


Substituting the values of G(s) and H(s) and simplifying , 


-8\ (44 S\(p4 & ros ) 
a. 0.154 (10 (1+ s\(1+ s) + Ont + 0.07 
nj | 
a am ane a! 
(1+ $)(1+ 35) + 0.1 (\+ sr) 
n 
the poles of }; determine whether or not the system is stable and these can be 


found by setting the denominator equal to zero and finding the roots of the resulting 
equation: 


(16) 


2 Ss _S_. oe = 
: (1+ 5) (1+ i) + 0.1(1+ sor) 0 (17) 


or, 


s 
aie ( ij — 


using the root-locus technique and an ESIAC potential-plane analog, we obtain the 
root-locus plot shown in Figure 6. 


=—| (18) 


Figure 6 shows that the poles of - are located in the stable region of the log s 
plane and their values are: 


P| = —35 

Po = —0.076 
P, = 1.52120 
P, = 1.5£-120 


The damping factor for the complex poles is found from the damping factor scale 
andis 0.5. 


on 
The zeros of fy, can be found by applying the above technique to the numerator 
of Equation (16). The equation to be solved is: 


6.5 (107) (1+ 535) 
(1+ $)(i+ $5) 


= = (19) 
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This is the same as Equation (18) except for the multiplying constant so the root- 
locus plot shown in Figure 6 can also be used to find the zeros. 
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FIGURE 6 LOG s PLANE ROOT LOCUS PLOT OF REACTOR CONTROL SYSTEM 
WITH ZERO TEMPERATURE COEFFICIENT 


As far as stability is concerned, there is no restriction on the location of the zeros 
and the values are seen to be: 


Z; = —0.07, Zo, 23, Za, AT 08) 


Using these pole and zero values, the complete transfer function with K;, =O1is: 


Ue. 
n 


(20) 


Ss ; Ss 
——— } (|\—- — i+ — 
1.52120) \ 1.52—120 c 


Notice that ng does not appear in the transfer function which means that the system 
behavior does not depend on the steady state power level. However, this may not 
be true when the temperature coefficient is other than zero. 


Let us now consider what happens to the control system if the tempe rature 
coefficient Kreis not zero. The complete reactor transfer function is given in 
Equation (13) and when it is incorporated in the control system, the overall closed 
loop transfer function is: 
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i 2 ()4 S\ (14 3 | $\(\4 8 _ 
Kas? (1+) (1+ + njKickrKa (I+ 3) (1+ G55) + 0.1 I+ ots | 
eat a 1 a PE BA oar Tee (21) 
2 Ss —s 
s? (1+ 3) (1+ ge) + nikicKrs (I+5 I+ oor + 0.1 I+ 37 


The intermediate steps are shown in the appendix. Notice that when Kj, =O, Equa- 
tion (12) reduces to Equation (16), 


The poles of this system are found by setting the denominator of Equation (2l) 
equal to zero and solving for the values of s that satisfy the resulting equation. 


nj Ke 1408 I+> =: 


—l = : (22) 


6.1 i——= |— : 
1.5Z120/ \ 1.52—120 


Here again, the easiest way to solve the above equation is to use the root-locus 
technique and an ESIAC potential plane analog. 


Figure 7 shows the root locus plot of Equation (22)as obtained with an ESLAC 
with the product njKi¢e as the varying parameter. Notice that the poles do not enter 
the unstable region as Kie is varied. However, one of the poles does approach the 
origin and for any reasonable value of Kt, can be considered to be at the origin. 
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FIGURE 7. LOG s PLANE ROOT LOCUS PLOT OF REACTOR CONTROL SYSTEM 
WITH NEGATIVE TEMPERATURE COEFFICIENT 
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The reference nj; should be approximately equal to nN) which is on the order of 
108 neutrons per cubic centimeter. This means that when Kj, is 10° or greater, the 
poles of the system will be at: 


$s; = .@) 
So = —2 
$3 —00 
which gives: 
n ZEROS 


Hoos ee (23) 
= (1+ 5) 


Although not shown here, the zeros are found in the same manner by setting the 
numerator of Equation (21) equal to zero and solving for the roots. The result gives: 


nL (- sha) (\—s37°75) 


"alteg) 


The fact that one of the poles of aq is at the origin is usually undesirable so the 
system should be modified. Notice that the poles of the system are identical to the 
poles of the control rod servo. This suggests that the system might be improved 
by changing the characteristics of the servo, as for example, by adding tachometer 
feedback. This modification is shown in Figure 5 by the dotted feedback loop around 
the control rod servo. The question immediately arises as to what value of tacho- 
meter constant K, should be used? 


(24) 


The control rod servo transfer function < becomes: 
2 


D _ | 
ee aT (25) 


s K, +1+8(I+3) 


which can be solved with a root locus plot. However, this is only a quadratic so 
using the value of K, that gives critical damping, the control rod servo transfer 
function becomes: (K; =1.83) 


cen eas ° 


This new transfer function can now be used in place of the original control rod ser- 
vo transfer function in the deviations of Equations (I6) and (20. The new pole loca 
The results are shown in Figures 8 and 9. Notice in | Figure 9 that when K,, differs 

from zero the poles of the system do not enter the unstable region or approach the 

origin. 
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FIGURE 8 LOG s PLANE ROOT LOCUS PLOT OF REACTOR CONTROL SYSTEM WITH 
ZERO TEMPERATURE COEFFICIENT AND MODIFIED CONTROL ROD SERVO 
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FIGURE 9. LOG s PLANE ROOT LOCUS PLOT OF REACTOR CONTROL SYSTEM WITH 
NEGATIVE TEMPERATURE COEFFICIENT AND MODIFIED CONTROL ROD SERVO 
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SUMMARY 


A potential plane analog of the log s plane can be used to great advantage for 
investigating problems involving reactor transfer functions and control systems, 
These uses include four fundamental operations: 


(1) Root Locus Plots 
(2) Bode Plots 

(3) Residue Evaluation | 
(4) Polynomial Factorization 


Although not illustrated in the above example, these operations can also be applied 
to problems involving: 


(a) Dead Time 

(b) Sampled-Data System 
(c) Linear Systems 

(d) Non-Linear Systems 


and others. 
The above example has illustrated the use of the potential plane analog for: 


root locus plots, (Figures 6, 7, 8, 9); Bode plots, (Figures 2, 3) and factoring 
polynomials, (Equations 18, 19, 22). 


REFERENCES 
(1) Morgan, M. L. "Algebraic Function Calculations Using Potential Analog Pairs", 
Procedure of the IRE, January 1961. 


(2) Looney, J. C. "The ESIAC petencalaniane Analog pe oa Automatic 
Control, Vol. 14, No. 5, p. 45, May 1961. 


(3) Morgan, M. L. and J. C. Looney, "Design of the ESIAC Algebraic Computer", 


IRE Transaction on Electronic Computers, Vol. EC~10, No. 3, September 1961. 


(4) Looney, J. C. “Sampled-Data System Design with a Potential Analog of the 
Log s Plane", Automatic Control, Vol. 15, No. 1, January, 1962. 


(5) Schultz, M. A. "Control of Nuclear Reactors and Power Plants", McGraw- 
Hill, 1955. 


EB-36 2/63 


«An 


APPENDIX 


DERIVATION OF REACTOR TRANSFER FUNCTION 
(Reference 5) 


6 
Since > Bi = B Equations (5) and (6) are combined to get: 
i=] 


dn _ nOk BS] Bi _ IG 
“dt. oS [7 


dt T = dt 

(Al) 
dn _ nAk _ S dC 
dt 7 & dt 


Now let n and Cj be separated into two parts; the steady-state values np and Cj,, and 
the small variations An and ACj. 


dn dAk _ nAk , Andk &dAC; 
oe a rn 2, dt ne) 
, dno dCio 
Since a. oF at @) 
and ene is small compared to a 
dC; dAC; No Bi An Bi 
also ao = Ta = a + —— — dj Cio — A, AC; 
. ; dAC; . Np Bi: (A3) 
which is at = Tr A; AC; 
nx Ri 
since dCi. = O = roe — i; Cio 
Now take the Laplace Transform of (A2) and (A3) to get: 
No Ak(s) 6 4 
sAn(s) = —— —§ >, ACi(s) (A4) 
i=l 
iAn(s 
sAC;(s) = = — d; AC;,(s) (AS) 


where initial conditions are equal to zero. Rearranging(A§5)gives: 


Bi An(s) 


Acie) = Ts4a) 


— dj; AC((s) 
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= Ge 


which can be substituted into (A4) and rearranged to get Equation (7). 


An(s) _ NG | 
Ak(s) = T | 
$ 


= B; 
+2 Ts | 


DERIVATION OF CLOSED-LOOP TRANSFER FUNCTION 


Nn = Not+An = not eG(s) 


e = nj —nH(s) 


Nn = no+ G(s) [ni — nH(s)] 
~= No+ nj G(s) — n G(s) H(s) 


; (1+-6(6) (3) atnnclal 


lle 


 (146()H() = $2 +6(s) & 1 +665) 


because nj 2 n and njp >> Anj. 


an. 1+ G(s) 
ni | + G(s) H(s) 
Referring to Figure 5: 
H(s) = Kg 


0 (6) Be 


(7) 


(15) 


EB-36 2/63 


SAG 


For Kie #0: 
en 2 noKR (I+ aor) 
Ak — : 
s(1+ +35 2 No Kr Kite ( = sir} 
Then: 
| Kr Kg Kr ( + at.) 
G(s) : 307 


Tika s2(14% cl) (ee 
No Kgs (1+ 5) (\+ Se} + nj Kic KaKg (1 + 5 \ (+ 5355) 
and the closed-loop transfer function is: 
2 Ss sits : = ioe! 
Kas (i+ $) (i+ <.) + nj Kte Kr Kas (I+ $) (1+ tots) + K, Kg Kr (1+ sia! 


n; 
s* (1+ 3) (i+ = + nj Kic Kes (1+ $)(1+ 585) + K, Ky Kp (I+ 557) 


which is Equation (2!) when the values of the constants are included. 


If Ky = O, the center terms in the numerator and denominator of Equation (2l) 
drop out and the result is Equation (16). 
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